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ABSTRACT 

In the course of a comprehensive mm/submm survey of massive star-forming 
regions, the vicinities of a sample of 47 luminous IRAS sources were closely in- 
vestigated with SCUBA and IRAM bolometers in order to search for massive 
protostellar candidates. A particularly interesting object has been found in the 
surroundings of the bright FIR source IRAS 07029-1215. Follow-up line obser- 
vations show that the object is cold, that it has a massive envelope, and that it 
is associated with an energetic molecular outflow. No infrared point source has 
been detected at its position. Therefore, it is a very good candidate for a member 
of the long searched-for group of massive protostars. 

Subject headings: stars: formation, ISM: jets and outflows, submillimeter 

^Max-Planck-Institut fiir Radioastronomic, Auf dcm Hiigel 69, D-53121 Bonn, Germany; forbrich@mpifr- 
bonn.mpg.de 

^Max-Planck-Institut fiir extraterrestrische Physik, Giessenbachstrafie, D-85748 Garching, Germany 



- 2 - 



1. Introduction 

The study of the earhest stages of massive star formation (MSF) is one of the most 
exciting subjects of present research in this field. Whether the dominant formation process 
for massive stars is disk accretion (see e.g. Jijina & Adams 1996) or coalescence (Bonnell, 
Bate, & Zinnecker 1998) is still an open question. The presence of outflows in regions of 
massive star formation favours, however, the first scenario. Similarly to previous surveys (e.g. 
Beuther et al. 2002, Mueller et al. 2002 or Hunter et al. 2000), we have searched for MSF 
candidates located near bright IRAS sources using SCUBA and IRAM bolometers (Posselt 
2003 and Klein et al. 2003, in prep.). Our focus was on objects unrelated to the IRAS 
sources themselves, looking for the earliest evolutionary stages not yet seen in the near- to 
mid-infrared wavelength regime. Even if we do not yet know what the observational features 
of massive protostars are (Evans et al. 2002), their invisibility at NIR/MIR wavelengths is a 
defining criterion. Near IRAS 07029-1215, an object with an IRAS luminosity of L{IRAS) = 
1700 Lq (Henning et al. 1992) at its distance old — 1 kpc, we discovered a deeply embedded 
object (UYSO 1, for Unidentified Young Stellar Object i), powering a high- velocity bipolar 
CO outflow. The IRAS source is located in the bright H ll region S 297, illuminated by the 
BlII/III star HD 53623 (Houk & Smith-Moore 1988). UYSO 1 is located very close to the 
edge of a dark cloud. 

2. Observations 

CO (J = 3 — > 2) line mapping of UYSO 1 was carried out with the James Clerk Maxwell 
Telescope (JCMT) ^ on Mauna Kea, Hawaii, in October 1999. The facility receiver B3 
(315-373 GHz) was used as frontend, the Digital Autocorrelation Spectrometer (DAS) as 
backend. With a total bandwidth of 250 MHz, centered on V;s,. = -|-12 km s~^, the maps 
were obtained in position- switch mode with an off-position 10' to the east. Map sampling 
was at intervals of 10", the total on+off integration time was 10 min, the beam efficiency 
r) — 0.63, and the beam size was 14". Additionally, SCUBA continuum observations at 
A = 450 /im and A = 850 /xm were performed with a total integration time of 46 minutes. 
Here, beam sizes were 8" and 14", respectively. Based on six different field scans with 
different chop throws, noise was minimised using the Emerson II technique (Emerson 1995). 
In March 2003, additional CS(J = 2^1), CS(J = 5^4), SiO(J = 2^1) and especially 



■'.JCMT is operated by the Joint Astronomy Center on behalf of the Particle Physics and Astronomy 
Research Council of the United Kingdom, the Netherlands Organization for Scientific Research, and the 
National Research Council of Canada. 
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H2C0( J = 3o3 — 202) /{J — 322 — 221) observations were obtained at the IRAM 30m telescope 
using the BlOO, B230, AlOO and A230 receivers. The raster map samphng was the same 
as for the JCMT observations, though a smaller range of offset positions was observed. 
The center position is shifted southward by 6.4" when compared to the JCMT map. Data 
was obtained in position switch mode, integration time at each offset was 6 min. All line 
observations together with the observing parameters are summarized in Tab. 1. 

We used the CLASS software developed by the Grenoble Astrophysics Group for the 
line data reduction, including zero-order baseline subtraction, and GAIA for the SCUBA 
data. 

3. Results 

The POSS'' view in Fig. 1 shows the environment of IRAS 07029-1215. The over- 
laid SCUBA map at 850 /xm displays UYSO 1 as a compact emission peak at a;(Bi95o.o) = 
07^02"'51', (5(Bi95o.o) = -12°14'26", between S 297 and the dark cloud. In addition to this, 
two dusty filaments are visible at 850 /um. The eastern filament contains the IRAS point 
source position, the western one is located inside the dark cloud. 

UYSO 1 is still invisible at near- and mid-infrared wavelengths: No infrared point source 
in the A =2. 2. .20 jim range is visible in 2MASS and MSX images (see Fig. 2a,b), but some 
faint "reflected" light can be seen in the 2MASS Ks-hand image. It cannot be ruled out, 
however, that this is captured H2 line emission, e.g. tracing cloud surfaces and outflows. 
Neither the IRAS, nor the 2MASS, nor the MSX infrared survey point source catalogues 
resulted in detections at the UYSO 1 position. There is, however, an MSX point source 
inside the "reflected" light south of UYSO 1. 

The CO map in Fig. 2c is a close-up view of the dense dust core. The line wings, 
extending from -20 km s"^ at [0",-20"] to +40 km s'^ at [0",0"] (see Fig. 3), can be clearly 
distinguished from well-determined baselines, extending from —40 km to -|-50 km s~^, 
respectively. Lines are narrower in off-outflow positions. Fig. 4 shows a comparison of 
the CO (J = 3 — >^ 2) line profile with the CS(J = 5 — > 4) profile in the map center. 
CS( J = 5 ^ 4) is a high-density tracer with a critical density of Ucnt = 7 x 10^ cm~^. 



"^The Digitized Sky Surveys were produced at the Space Telescope Science Institute under U.S. Govern- 
ment grant NAGW-2166. The images of these surveys are based on photographic data obtained using the 
Oschin Schmidt Telescope on Palomar Mountain and the UK Schmidt Telescope. The plates were processed 
into the present compressed digital form with the permission of these institutions. 
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Contrary to what could be expected from the 850 //m map, the CO data do not show 
a spherically symmetric structure. Instead, the envelope seems to be compressed towards 
the Hll region east of it. The CO peak is roughly 10" off the SCUBA peak position. The 
line wing map displays a bipolar outflow in the NW-SE direction, its origin being UYSO 1 
at raster coordinates [0",-10"], taking the emission peak of the complete integrated line 
as definition. The outflow is not clearly visible in the CS data, possibly due to the average 
density being below the CS critical densities. Lines on corresponding positions have, however, 
corresponding shapes. SiO(J = 2 — > 1) could only be detected at the position [-10", 16.4"], 
in the outer regions of the redshifted outflow lobe. Possibly, this is due to a bow shock caused 
by the outflow coUiding with the local interstellar medium. The detection at this position is 
a 6 (T detection. With the spatial resolution of the present observations, we cannot decide 
whether the bipolar outflow is composed of only one or more components (cf. Beuther et al. 
2002b). 

The SCUBA 450 fxm map (see Fig. 2d) only shows a point source 6" to the north of 
UYSO 1, but in between the two outflow lobes. 

4. Analysis 

Prom the CO (J = 3 — > 2) observations, the object mass and outflow parameters can be 
deduced. The mass derived from an integration of column densities can be compared to the 
virial mass. Line wing analysis delivers physical parameters of the bipolar outflow, however 
still dependent on the unknown outflow inclination angle. The distance to UYSO 1, assumed 
to be about the same as the distance to IRAS 07029-1215, was kinematically determined to 
be 1 kpc by Wouterloot & Brand (1989). However, while the distance to the illuminating star 
of S 297, namely HD 53623, was spectroscopically determined to be 2.2 kpc ± 25% by Neckel, 
Klare, & Sarcander (1980), Hipparcos measured a parallax for HD 53623 of tt = 2.21 ± 1.01 
mas, corresponding to a distance of 4501i4q PC (ESA 1997). Assuming the true distance 
is near the upper limit of this interval, we adopt a value of 1 kpc for the following mass 
estimates, in accordance with detailed studies of CMa Rl by Shevchenko et al. (1999) and 
de Zeeuw, Hoogerwerf, & de Bruijne (1999). 

The CO data analysis is depicted in Pig. 2c, in direct comparison to 2MASS-ir5^ data 
and the 450 /im continuum. 



^ Atlas Image obtained as part of the Two Micron All Sky Survey (2MASS) , a joint project of the University 

of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded 
by the National Aeronautics and Space Administration and the National Science Foundation. 
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4.1. Temperature Determination 

While the main beam temperatures Tmb of the highly optically thick CO lines already 
confine the kinetic temperature to the 40-50 K range, another estimate was obtained from 
H2C0(J = 3o3 ^ 2o2)/(J = 322 ^ 22i) data, following Mangum k Wootten (1993). Fig. 5 
shows the two H2CO transitions for UYSO 1. The intensity ratio deduced from fitted Gaus- 
sians is 6, thus indicating a kinetic temperature of T=40±5 K for a relatively wide range of 
densities. 



4.2. Mciss Determination 

Assuming virial conditions, i.e. assuming that the observed line width is entirely due 
to motion in gravitational equilibrium, the total mass of the object can be estimated as 
described by Henning et al. (2000). Inside a sphere of 0.14 pc radius, corresponding to the 
40% contour level (deconvolved with the 14" beam) at 10 cr, and using the Gaussian fit result 
of Av — 3 km s~^, a virial mass of Myi^ — 180 M0 is obtained. 

Following Henning et al. (2000), we derive hydrogen column densities from the CO(J = 
3 — > 2) transition, assuming a relation of A'"(H2) (cm~^) = 3 x 10^° J Tmbdv. An under- 
lying assumption is a conversion factor of CO (J = 3^2) / CO (J = 1 — > 0) close to 
unity, as it is frequently observed in massive star-forming regions and explained by clumpy 
structure (Storzer et al. 2000). The gas mass of the cloud can then be determined as 
Mgs.s = 2mH/iA(A^(H2)), where /2 = 1.36 takes into account the interstellar mean molec- 
ular weight per H atom. A is the area occupied by the source. 

Within the same area and to the same confidence level as the virial analysis, this leads 
to a different result for the cloud mass. With an average column density of A^(H2) = 
3.2 X 10^^ cm"^ the resulting mass is Mgas = 40 Mq. 

Another estimate for the total mass has been calculated from the continuum obser- 
vations at A = 850 fim. Assuming a dust temperature of = 20 K and optically thin 
conditions, Mgas = F^, ■ (P /{Km{i') • Bi,{T^)) • [Mgas/Mdust]- Using an interpolated opacity of 
K = 1.97 cm~^/g (thin ice mantles, Ossenkopf & Henning 1994), the result is M = 15 Mq. 
The radius of the region taken into account for this calculation is 0.11 pc (at a distance of 
1 kpc), slightly smaller than the 0.14 pc used here. Additionally, the dust continuum emis- 
sion traces only the high-density central regions while the CO emission is collected along 
the line of sight throughout the low-density halo. Optically thick emission might further 
infiuence the estimate (see §4.4). These two points at least partly explain the discrepancy. 
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Thus, the mass estimates based on line and dust emission are compatible when keeping 
in mind their respective limitations and range between 15 and 40 Mq. The virial mass is 
much higher which is certainly caused by an overestimation of line width due to optical 
thickness and turbulent motion, in addition to the cloud not being in virial equilibrium. 

4.3. Outflow Properties 

The two line wings were analyzed to their respective 10% contour levels, corresponding 
to a 2 o" detection for the redshifted line wing and a 4 cr detection for the blucshifted line 
wing. The masses inside the two outflow lobes are again calculated by integration of H2 
column densities, resulting in Mbiue = 1-8 Mq and Mred = 3.6 Mq, the total outflow mass 
thus being Mout = 5.4 M©. 

The maximum projected velocities of the red and blue line wings are t'proj.r = 27.9 km s~^ 
and Vpj-o},h = 32.1 km s~^ respectively, the mean value being Wpj-oj = 30 km s^^. Together with 
the known distance and the outflow radial sizes, the dynamical timescale and the mass loss 
rate can be obtained. With caution, the latter can be empirically converted to luminosity, 
and spectral type as well as mass of an assumed central star. All values discussed hereafter 
are summarized in Tab. 2 for inclinations oi i — 10° i — 57.3° (the most probable value, cf. 
Bontemps, Ward- Thompson, & Andre 1996) and i — 80? 

When the apparent outflow velocity Vproj and the apparent maximum outflow extension 
have been inclination-corrected to Vout and -Rout, the dynamical timescale can be calculated 
from td{i) = -Rout(0/'^out(0- The resulting values range from id(80°) pa 500 years to id(10°) ~ 
1.5 X 10^ years. However, the timescale estimated this way is not a good age indicator. 
Only a massive object can uphold mass outflow rates M{i) — Mout{i)/td{'i) ranging from 
M(10°) ^ 4x 10"^ Mg/yr to M(80°) ^ Ix lO'^ Mo/yr. They can be empirically transformed 
into the luminosity of the central object according to Shepherd & Churchwell (1996), Henning 
et al. (2000) or, based on a more comprehensive sample, by Beuther ct al. (2002b). Tab. 3 
summarizes the resulting spectral types, masses and luminosities (using data by Lang 1991). 
However, these properties of the central object can be better determined using the SED (see 
§4.4). 

It may be interesting to consider the age of the neighbouring Hll region. Following 
Cernicharo et al. (1998), this age can be estimated (Whitworth et al. 1994). The size of S 297 
in the optical regime with a diameter of 7' (Sharpless catalogue) corresponds to a radius of 
only about 2.1 x 10^ AU, hence indicating that the nebula is rather young. Assuming a 
typical mean density of 10^ cm~^ and a Lyman continuum flux for a Bl star from Thompson 
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(1984), the resulting age is approximately 7 x 10^ yrs. 

Regardless of the uncertainties in these age estimates, the H ll region is hkely to be older 
than UYSO 1. Possibly, this is an example of induced star formation, although this would 
need kinematic confirmation. 

4.4. Spectral Energy Distribution 

While the two measured continuum fluxes at 450 /im and 850 /xm determine the cold dust 
emission, more data had to be complemented before the spectral energy distribution could 
be tentatively compared to modeling results. Given the aim of determining the luminosity 
and mass of the central object as well as having a comparison to the mass determination, 
upper flux limits for short wavelengths helped pinning down physical properties of the central 
object. As main support besides the JCMT data, IRAS fluxes were used, determined at the 
position of UYSO 1 within beam sizes from the IRAS atlas maps. UYSO 1 is unresolved 
also on IRAS HIRES maps. Additionally, data for the definitely brighter nearby MSX point 
source G225. 4582-02. 5939 and 2MASS detection limits set further constraints. 

Integrating the SED and using the fiux limits at infrared wavelengths as the actual 
fiux values leads to a luminosity of < 1900 Lq, assuming spherical symmetry. A modified 
blackbody accounting for the dust optical depth was fitted to the data: S,y = B,,{T){1 — 
exp(— rj,))Af2, with the beam size Afl and with dust opacities Hi, oc (see e.g. Hildebrand 
1983). A good fit is obtained for T = 45 K and (3 — 2, however fits with temperatures 
between 40 and 50 K are also fitting the data within the uncertainties. With the aim of 
checking the consistency of existing data, the continuum emission was then modeled using 
the radiative transfer code developed by Manske, Henning, &; Men'shchikov (1998). This is 
an accelerated version of the 2-D ray-tracing code developed for radiative transfer in disk 
configurations by Men'shchikov & Henning (1997). Due to insufficient parameter constraints, 
the code was used in its ID version only. 

The limiting outer radius of the model could be estimated from the CO data to be 

3.2 X 10^ AU, taking the same area as for the mass estimation. Assuming an exponential 
of the density distribution of 2.0 with a dust-to-gas ratio of 1:150, the model uses silicate 
optical data from Dorschner ct al. (1995), while the properties of the carbonaceous dust are 
from Jager et al. (1998) [1000 K data]. The ratio of silicates to carbon is 3:2. The inner 
model radius is determined by the dust sublimation temperature of 1500 K. 

Fig. 6 shows the results for an envelope mass of Mcnv = 30 Mq and an assumed single 
central object in the B2.5 class. The optical depth at A = 500 nm for the simulation runs 
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is T5oonm ~ 6400. Since the code becomes inaccurate for high optical depths, the computed 
SED is arbitrarily shown for r < 1000, corresponding to A > 10 /xm. At A = 850 /xm, the 
optical depth is rgso/im ~ 1- Thus, the mass estimate from the dust emission in §4.2 is rather 
a lower limit because of contribution from optically thick emission. 

Both the simulation and the blackbody fit indicate the equivalent of an early B star 
as the central object. This finding is compatible with the result of the empirical outflow 
analysis for i = 57.3°. The envelope mass in the simulation is compatible with the mass 
estimate from the CO(J = 3^2) data (40 Mq), especially when keeping in mind that the 
modeled gas density only reaches 10"^ cm~"^ at its outer radius, thus allowing for a slightly 
higher envelope mass when integrating further outwards. 

5. Conclusions 

The emerging picture of UYSO 1 drawn from both mass estimates and subsequent 
empirical relations as well as SED considerations thus is that of an early B star surrounded by 
an envelope of about 30-40 M0. The object seems to be in a particularly early evolutionary 
state, but already drives a bipolar outflow with a total mass of Moutfiow = 5.4 Mq. To 
our knowledge, there are only very few comparable sources that are undetected at infrared 
wavelengths (see e.g. Hunter et al. 1998, Bernard, Dobashi, & Momose 1999 and Fontani et 
al. 2003). With additional high-resolution observations, we hope to be able to specify these 
conclusions and learn more about this evolutionary state. Possibly, parameters of a putative 
disk around UYSO 1 could be derived then. 

The authors are grateful to Remo Tilanus for his support during the JCMT observing 
run, as well as to R. Mauersberger for carrying out the DDT observations at the IRAM 30m 
telescope in service mode. 

This research has made use of the SIMBAD database, operated at CDS, Strasbourg, 
Prance. 
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Table 1: Heterodyne observations of UYSO 1 



Line 




Frequency 


Telescope 


Date 


beam size 




C0(J = 3 


-2) 


345.79 GHz 


JCMT 


Oct 99 


14" 


0.63 


CS(J = 2- 


-1) 


97.98 GHz 


IRAM 30m 


Mar 03 


25" 


0.78 


CS(J = 5- 


-4) 


244.94 GHz 


IRAM 30m 


Mar 03 


10" 


0.50 


SiO(J = 2 




1 86.84 GHz 


IRAM 30m 


Mar 03 


29" 


0.78 


H2C0(J = 


3- 


2)^ 218.34 GHz 


IRAM 30m 


Mar 03 


12" 


0.56 



-^HaCOCJ = 3o3 2o2) and (J = 822 ^ 221) 



Table 2: UYSO 1 outflow properties 



i 


Vout 


R-out 




td 


M 




10° 


30.5 km s-^ 


1.4 X 10^^ 


km 


14600 yrs 


3.7 X 10-^ 


Me/yr 


57?3 


55.5 km 


2.9 X 10^2 


km 


1700 yrs 


3.3 X 10-=^ 


Mo/yr 


80° 


172.7 km s-^ 


2.4 X 10^2 


km 


450 yrs 


1.2 X 10-2 


Mo/yr 



Table 3: Outflow-inferred properties of assumed central star 



i 


L 


spec, type 


Mgtar 


10° 


10^ Lo 


B4 


6.5 Mq 


57?3 


10^ Lq 


Bl 


13 Mq 


80° 


10^ Lq 


09 


19 Mq 
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Fig. 1. — The surroundings of UYSO 1: POSS-I, overlaid with 850 //m continuum from the 
JCMT. The box indicates the FOV of Fig. 2. 
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7^02"'52^ 50^ 30 -30 
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Fig. 2.— Close-up views of UYSO 1 in a) MSX-A, b) 2MASS-X5, c) CO(J = 3 -> 2) 
and d) 450 /xm continuum (both JCMT). The box in the MSX data shows the FOV of the 
other pictures. The integrated CO hne data is shown in 10% contours, the two hne wings 
are shown down to 50% intensity. Spectra positions are indicated as well as the beam sizes 
for CO and the 450 //m continuum. UYSO 1 is located just below the center position at 
[0",-10"], the position of the CO peak. 
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Fig. 3. — C0( J = 3—^2) spectra on the inner grid. The hne maxima (aU between 40 K and 
50 K) were cut to allow for a better depiction of line wings. 
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Fig. 4. — Comparison of the C0( J = 3^2) and CS(J = 5^4) lines at the respective peak 
positions: the CO hne at [0",0"] and an addition of four CS spectra, amphfied ten times for 
depiction, at [0",-6.4"]. Positions are in arcseconds relative to the center position of the CO 
grid. 
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Fig. 5. — Temperature determination with H2CO lines (see text). Gaussian fit results are 
shown in grey. 
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Fig. 6. — Modeling results from the Manske code. The straight line shows the total flux 
(displayed for tx < 1000), boxes indicate beam-convolved fluxes. The dashed line corresponds 
to a modified blackbody at T=45 K (with /? = 2) 



